Effect of viscous dissipation and induced magnetic field on an unsteady mixed convective stagnation point flow of a nonhomogenous nanofluid

In the study, we investigate the numerical investigation of variable viscous dissipation and source of heat or sink in mixed convective stagnation point flow the unsteady non-homogeneous nanofluid under the induced magnetic parameter. Considering similarity conversions, the governing of fundamental boundary of layer non-linear PDEs are transformed to equations of the non-linear differential type that, under appropriate boundary conditions, are numerically solved, and the MATLAB function bvp4c is considered to solve the resulting system. The obtained results are calculated numerically for non-dimensional velocity, temperature, and volume fraction and displayed graphically. Further, numbers of Nusselt and Sherwood and local Skin of friction have been produced and displayed by graphs. A comparison with previous results obtained neglecting the new parameters has been made to show the impact of new external parametes on the phenomneon. The obtained findings agree with those introduced by others if the magnetic field and viscous dissipation are neglected. The results obtained have an important applications in diverse field as chemical engineering, agriculture, medical science, and industries.

Moreover, the colloidal mixture of particles nanosized in a fluid base enhances transferring heat suited features in real applications.Dharmalingam et al. 1 presented a mathematical and experimental study about transferring heat through nanofluids, as well as the chemical and physical features, and analyzed the opportunities and challenges of nanofluids in further studies.Furthermore, they reviewed nanofluid, the latest research, its important heat transfer mechanisms, and development in nanofluids.Even though the nanofluid usage in the experimental application is critical due to forming sedimentation and clogging in the flow path, many of the mechanisms of transferring heat are discussed.Certainly, the use of liquids in nanofluids requires an inclusive understanding of their properties, such as thermal conductivity, specific heat, viscosity, etc.To expand nanofluid applications, we need further studies on the contributing factors for heat transfer and friction.The presence of suspended nanoparticles that have high thermal conductivity increases the higher thermal capacity of the liquid.
The flow-through boundary layer on the surface that moves continually of the 2D domain was at first modelled by Ref. 19 .The heat transfer and mass transfer convective process takes an important place because of the effects of buoyancy because of different concentration and temperature profiles, in turn.An important role has been played to address the phenomena transport, the thermal and mass diffusion taking place with the simultaneous buoyancy forces action.Mohamed et al. 20 examined the finite element method technique (FEM) of a flow hydromagnetic and transferring heat of the fluid generation of heat over the embedded surface within the non-Darcian porousity considering chemical reactions.In Ref. 21 , the authors pointed out the chemical reaction and thermal radiations effect on transferring heat and mass in MHD flow of micropolar over a moving vertical porous plate within a porous medium that has a generation of heat feature.Effect of MHD and thermal radiation free convective flow of a polarized fluid through a porousity was discussed (see 22 ).Analytical solution of the chemical reactions and thermal radiations impacts on MHD convection unsteady in a porous medium under heat source/sink were explored Abdel-Nasser et al. 23 .Mohamed et al. 24 pointed out the unsteady magnetohydrodynamic convection boundary layer flow double-diffusive past a radiate vertical surface hot within a porous medium in considering chemical reactions and sink of heat.In Mohamed et al. 25 , the authors studied the chemical reactions, thermophoresis, variable viscosity, and heat source/sink impacts of on the non-Darcian mixed convective heat, as well as the flow of mass transfer.Carreau nanofluid radially stretching surface under effects of the magnetic field and the partial slip on an unsteady axisymmetric flow were discussed by Azama et al. 26 .Recently, the authors of Khashi'ie et al. 27 highlighted the heat transfer and mixed convective flow of the dual stratified micropolar fluid produced using the permeable shrinking/ stretching sheet.Kumar and Srinivas 28 discussed the Eyring Powell nanofluid flow unsteady hydromagnetic over stretching sheet inclined permeable under Joule heating and thermal radiation.
Nanofluid with magnetic field represents a colloidal suspension stable containing the base the solution, magnetic particles single-domain and surfactant.Because the motion of Brownian of the particles and the repulsion of the surfactant fosters the magnetic fluid and the resistance of magnetic particles to gravity sedimentation, is stable and depends on the size of magnetic nanoparticles and the physical properties of the based fluid (see Puga et al. 29 ).Thus, thermal management and heat transfer have garnered widespread attention, because of their particular electronics devices, industrial applications, and thermal.The magnetic Fe3O4@CNT nanocomposites, the viscosity, thermal conductivity, and capacity of specific heat were determined by Shi et al. 30 .The adjustability, magnetic responsiveness, and recyclability of the transfer of heat demonstrated experimentally, and the magnetic nanofluid was prepared by the properties thermophysical of the nanofluid and hydrothermal method based on the temperature were explained 31 .The transpiration andiscous dissipation effect on the generation of entropy inflow of hybrid nan-ofluid over a radially non-linear stretching, the disk is investigated by Farooq et al. 32 .Ahmed et al. 33 discussed the detail of the dual solution in a boundary layer flow of a fluid power-law over a moving plate permeable flat with viscous dissipation, heat generation/absorption, and thermal radiation.The combined electric field and magnetohydrodynamic impact on Maxwell's nanofluid flow an unsteady over a stretching surface under considering the thermal radiation and heat effect variable explained by Khan 34 .The efficiency of heat exchange of magnetic nanoparticles suspension, magnetic nanofluid, in natural convection can be adjusted by changing the direction, form, and strength of magnetic field due to its properties of excellent magnetic (see 35,36 ).However, the magnetic nano-particles thermal conductivity less than most nano-particles as CuO and Al2O3, limiting their wide applications in particularly in heat exchangers and the energy industry (see 37,38 ).Hence, nano-composites of magnetic with a combination of magnetic property and higher thermal conductivity are an excellent choice for increasing the exchange of heat efficiency of natural convection more than other varies materials strongly.Reference 39 discussed in detail the heat transfer convective magnetocontrollable of nanofluid via a straight tube.
Nanda et al. 40 discussed the enhanced heat transmission in methanol-based AA7072/AA7075 tangent hyperbolic hybrid nanofluid flow along a nonlinear expandable surface.Kumar et al. 28 explained: A nanotechnology application considering the effect of electromagnetic induction on the heat transmission in engine oil-based hybrid nano and ferro-fluids.Sandeep et al. 41 discussed the low and heat transfer in radiative MHD dusty-hybrid ferrofluids.Wang et al. 42 discussed the effects of nanoparticle aggregation and radiation on the flow of nanofluid between the gap of a disk and cone.Umavathi et al. 43 investigated magnetohydrodynamic squeezing Casson nanofluid flow between parallel convectively heated disks.Kumar et al. 28 studied heat transfer analysis in threedimensional unsteady magnetic fluid flow of water-based ternary hybrid nanofluid conveying three various shaped nanoparticles: A comparative study.Oke et al. 44 exploration of the effects of Coriolis force and thermal radiation on water-based hybrid nanofluid flow over an exponentially stretching plate.Gowda et al. 45 pointed out three-dimensional coupled flow and heat transfer in non-newtonian magnetic nanofluid: An application of

Mathematical analysis
This paper is a two-dimensional steady of coupled heat transfer by mixing the convective flow of incompressible viscous laminar fluid.Authors argue that the fluid is an unsteadily convective mixed stagnation point flow of nonhomogenous nanofluid and its various properties because of the restricted viscosity, density, and temperature.The variation density, and the buoyancy effect are considered in the equation of momentum (Boussineq's approximation), and species of concentration far from the wall is C ∞ , is infinite similarly small and hence neglecting the Sort and Dufour numbers.Consider the x-axis in the plate direction, and the y-axis is normal on it.The external flow in the direction parallel to the plate inclined and has a uniform velocity U ∞ .The surface temperature becomes uniform T w that is more than the ambient temperature T ∞ .The species concentration at the surface is uniform at wall C w , that is more than the ambient spices concentration C ∞ .The source/sink term of the chemical reaction is involved in the concentration equation.
The derivation of magnetic term is introduced as follows: The law of Ohm is written as (see 48,49 ).
where E represents the intensity electric field, j denotes the density of current, B = (0, B • , 0) , E = 0 , w e τ e = m and j y = 0 .Figure 1a shows the magnetic field in the coordinates system.We have Substitute from Eqs. ( 2) and (3) into Eq.(1); yields to Equation (4) implies that Hence, in the x-momentum equation, we add the force generated by the magnetic field, F B (1) www.nature.com/scientificreports/From Eqs. ( 6) and ( 7) yields w = 0 since there is no flow in z direction, Since τ e << and this is a result of a small magnetic field, the medium itself makes the electric conductivity small and m = 0 , therefore, Eq. ( 9) takes the form:

Mathematical formulation
We take into account the unsteady 2D fused convective boundary layer flow of an electrically conducting viscous fluid near the recession point on a vertical plate.The fluid electrical conductivity is due to nanoparticles suspended in the fluid.This paper studies the transfer of heat and fluid flow near the wall, where point of stagnation occurs because of the flow divergence.Figure 1b shows configuring the problem.The governing equations are Continuity equation: The equation of momentum: Energy equation: Concentration equation: In Eqs. ( 11)-( 14), u and v denote the velocities of the component along the axis x and y direction, in order.ρ f denotes the base fluid of density, T represents temperature of the nanofluid.σ and α represent the electrical of fluid and thermal conductivity, in order.β denotes the coefficients of thermal expansion, the thermophortic diffusion coefficient denoted by D T , D B suggests the coefficients of Brownian diffusiont.µ f represents the vis- cosity of the nanofluid, ρ p describes the nanoparticles density.φ denotes the nanoparticle concentration and g denotes gravity-caused acceleration.
The generalized Bernulli's equation, in free stream, in Eq. ( 12) is Substituting Eq. ( 15) into Eq.( 12) yields: Changing the nonlinear differential equations system into a set of ordinary differential equations system, the dimensionless variables become 50,51 , ( 7) where the ψ is the function of stream defined as u = ∂ψ ∂y , v = − ∂ψ ∂x .Thus, Eq. ( 11) in the model is fulfilled.The velocity components u and v have the following forms:

Boundary settings
The initial and boundary conditions are The ambient fluid velocity far away from the considered wall and determines the stagnation point flow strength, ct < 1 when a > 0 .Here, a and c are the constants that have dimension 1/time.The plate surface has stretching /shrinking along the velocity u w (x, t) = bx 1−ct where b < 0 represents the shrinking of velocity and its stretching velocity when b > 0 and c has dimension 1/time.The nano- fluid is non-homogeneous with average physical features of basic fluid and nanoparticles.At the plate x = 0 , the temperature T and the nanoparticle fraction φ fixed values T w and φ w , in order.The ambient values, attained as y becomes infinite of T and φ given by T ∞ and φ ∞ in order.
-The temperature of the wall takes T ∞ and T 0 denote the temperature of ambient fluid and the characteristic temperature, respectively.
-Further, we assume that and φ w (x, t) = φ ∞ + φ 0 x (1−ct) 2 where the term V * w denotes the uniform surface of mass flux.It is worth noting that V * w > 0 refers to injection wheres V * w < 0 means suction and it is an impermeable surface when V * w = 0 .A rigid vertical wall is located at x = 0

Numerical method
Employing the similarity transformation (17), Eqs. ( 13), ( 14) and ( 16) decline to these nonlinear ordinary differential equations: In this section, the MATLAB function bvp4c is applied to solve our model with boundary conditions.
The boundary settings become: denotes the nanofluid ratio of buoyancy, P r = µ f α f suggests the Pradentle number, Sc = ν f D f represents the nanofluid Schmidt number.
is the parameter of thermophoresis.
The utilized physical quantities denote the coefficient of the friction of skin C f , Nusselt and the Sherwood numbersNu x and Sh x that take the form as follows: In which the mass flux at surface q m and wall heat flux q w take the form Substitution Eqs. ( 17) and ( 18) into Eqs.( 25) and ( 26), we obtain: where Re x = U ∞ x ν f suggests the Reynolds number.

Results and discussion
The numerical solutions of Eqs. ( 21)-( 23) within the boundary settings (24) are displayed in Figs. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 and 15.The findings of the coefficient of local Nusselt number and skin friction are comparable with those of Ramachandran 52 , Lok et al. 53 , and Kashani et al. 54 for = 1, and various values of Pr, neglecting the parameters Scn, Nb, Nt, Nr, A, M, Ec, Vw, and ε .We noticed that the comparison shows a perfect agreement, and our approach gives identical results with the results obtained by others, so we concluded that the current results obtained are appropriate.Tables 1 and 2 confirm the validity of the current results for the Nusselt number and coefficient of skin friction.Firstly, the velocity f ′ starts from its minimum value (i.e., f ′ = b a = ε) at the wall (i.e., η = 0)and increases approaches to the unity as η → ∞ (far from the wall).The volume fraction φ and profiles of temperature θ begin with the unity at the wall and decrease approaching zero as η → ∞.This physically indicates the increasing of velocity but vanishing temperature and volume fraction when the fluid is far from the wall.The result agrees with the experimental findings.
Figure 2 represents the thermophoresis impact Nt on the profiles of velocity f ′ , volume fraction φ and tem- perature θ for the values Scn = 0.22 , Nb = 0.2 , M = 5 , Ec = 1 and A = Nr = Vw = 0.5 A significant impact is noticed on the boundary layer regrading various values of Prandtl number Pr when it takes the values of 0.72 and 9.In this case the velocity f ′ declines with the increased values of thermophoresis Nt, but the temperature θ and the volume fraction φ profiles rise.On the other hand, we show that the velocity declines with the increased values of Pr but volume fraction and the temperature decline.Similarly, Fig. 3 displays the impact of Brownian motion Nb for the velocity f ′ , the volume fraction φ , and the temperature θ profiles.No significant impact on the temperature, velocity and volume fraction when the Pr takes the values 6.2 and 9 respectively.Only when the Pr takes the values 0.72 and 9 then the temperature and velocity increase with the increased values of Pr and Nb, but the volume fraction decreases.This indicates the the same impact of Pr and Nb on all physical quantities.Figure 4 displays the impact of parameter of unsteadiness A on the velocity in x direction when = −2.2 .In this figure, we notice that the velocity is close to zero at the wall.When = −2 decreases, but when = 2 increases.On the other hand, the temperature and the volume fraction profiles decline concerning the negative and positive values of .Physically, this means to the positive impact for the positive value of but negative for the negative value of and has applications in the related phenomena.Figure 5 demonstrates the velocity profile behavior in the presence of ε parameter (the stretching ratio of the sheet) and A (unsteadiness parameter) for two cases: one for a steady A = 0 and other for A = 5 .The ε < 1, the velocity increases rapidly to reach the level of η .After that, it is identical with increase of η .In the same manner, when ε > 1, the velocity decreases dramatically and then becomes uniform with increase of η .Furthermore, the impact of ε parameter and A on the temperature unsteadiness parameter can be shown in this figure.One can demonstrate that the volume fraction and the temperature decline with increasing ε as well as the presence of parameter of unsteadiness A.
Figure 6 illustrates the parameter of magnetic M effect on the velocity, temperature and volume fraction profiles.The temperature only changes near the wall and then it stays uniform through the fluid as it increases.The thickness boundary layer is not under the impact of the magnetic field, due to the volume fraction profile.Therefore, the magnetic field affects the regions that have high concentration of the nanoparticles.The changes near the wall lead to the increment the transfer of heat, which can result in the change in the nanoparticles.Therefore, the temperature changes.Figure 7a plots that higher the parameter of the magnetic parameter is, the lower the velocity away from the wall becomes.Physically, we can conclude that when applying the magnetic field inside the boundary layer, it produces a type of resistance force called Lorentz body force that reflects the flow and slows the movement of the fluid.Figure 7 depicts the impact of Eckert number Ec on the temperature, the velocity, and the heat transfer.Ec statistically affects the fluid temperature because of improving Ec number, the improved viscous dissipation impact on the boundary layer of thermal.The Ec has a strong effect on velocity that we can deduce that the velocity along the surface increases and then becomes steady as η increases.Furthermore, we can see that there is a small change in the boundary layer thickness on the transfer of heat near the wall and then the boundary layer declines to become uniform (no change).
Figure 8 plots the thermophoresis effect Nt and the Eckert number Ec on the friction of skinCf Nu x also, the Sherwood number Sh x .Nt and Ec rise, the skin friction along with Nusselt number decline.In contrast, the Sherwood number rises with rising Nt but declines with the large values of Ec.
Figure 9 displays the impact of Eckert number Ec and Vw on the friction of skin Cf Nu x and Sherwood num- ber Sh x .Cf Nu x rises when rises when Ec and Vw rise, while the Nusselt number and Sherwood number decline with the increasing of Ec and Vw.Figures 10 and 11 display 3D of parameter of magnetic M variation and η, as well, variation of thermophoresis Nt and η on the nondimensional velocity f ′ the surface, respectively.
It can be seen in the two cases that the magnetic M and thermophoresis Nt enhance the velocity near the wall.Then, it becomes a constant away through the surface as M, Nt and η increase.Figure 12  various values of local Nusselt, friction of Skin, also, Sherwood number with magnetic M and mixed convection parameter regarding a fixed value of Ekert number Ec at 1.The friction of Skin rises with increasing the values of the parameter of mixed convection .On the other hand, there is a significant reduction in the Nusselt number and Sherwood number or heat transfer when the parameter of mixed convection rises with the parameter of magnetic field.Figure 13 displays the influence of wall transpiration parameterVw and Eckert number Ec in the (a) Skin friction, (b) local Nusselt number, and (c) Sherwood number.Skin friction decreases with rising the values of Vw and increase Ec, but the local Nusselt number and Sherwood number rises when the values of Vw and Ec enhance.Furthermore, it is seen that Fig. 14 exhibits the effect of unsteadiness parameter A and parameter of wall transpiration Vw on the Nusselt number −θ ′ (0) .Increasing A andVw rises the local Nusselt number.Figure (15) reveals that the positive values of the parameter of fused convection and parameter of unsteadiness A cause a reduction in the Nusselt number −θ ′ (0) .Nu reduces with rising A, but declines with increasing the fused convection parameter .
The obtained findings agree with those introduced by others if the magnetic field and viscous dissipation are neglected.The results obtained have important applications in diverse field as a chemical engineering, agriculture, medical science, and industries.Finally, if the viscosity is neglected, the obtained results in this study reduce the results obtained by Hamad et al. 9 .Furthermore, if the magnetic field induced vanishes and considering the flow of hybrid nanofluid over the non-linear radially enlarging disk.Those obtained findings deduce those of Farook et al. 32

Conclusion
We examined how the magnetic field and the viscous dissipation parameters affected the 2D unsteady mixed convective boundary layer flow.We found some vital point in the obtained results, such as: -Temperature surface θ as well as volume fraction φ increase when thermophorsis Nt increases.
-The dissipation of viscous Ec strongly impacts temperature θ by increasing its profile.www.nature.com/scientificreports/-The parameter of magneticM strongly impact on the temperature θ and f ′ near the wall.
-Eckert number Ec and thermophorsis Nt reduce the skin of friction C f and Nusselt number Nu, but raise the Sherwoord number Shx.
-The paper provides a system for achieving a superior performance of viscous dissipation variable and a heat source or sink on the unsteady mixed convective stagnation point flow of a non-homogeneous nanofluid under the induced magnetic field, and illustrates a system of improving the transfer of heat in a heat exchanger by controlling the magnetic field parameter distribution.-The temperature, velocity, volume fraction, Sherwood number,as well as Nusselt number are affected strongly by the magnetic field and viscous dissipation in an unsteady mixed convective stagnation point flow of the nonhomogenous nanofluid that indicated the applicability in diverse fields, such as medical engineering, nanotechnology, and petroleum extracting.
Our results are generalized, and the agreement of the results obtained by others was highlighted 9,32 with neglecting the new impacts.= f ′′ (0) Pr Ramachandran et al. 52 Lok et al. 53   = −θ ′ (0) Pr Ramachandran et al. 52 Lok et al. 53

Figure 1 .
Figure 1.Geometrical shapes of (a) the magnetic field through fluid (b) flow configuration with coordinate system.

Figure 6 .
Figure5demonstrates the velocity profile behavior in the presence of ε parameter (the stretching ratio of the sheet) and A (unsteadiness parameter) for two cases: one for a steady A = 0 and other for A = 5 .The ε < 1, the velocity increases rapidly to reach the level of η .After that, it is identical with increase of η .In the same manner, when ε > 1, the velocity decreases dramatically and then becomes uniform with increase of η .Furthermore, the impact of ε parameter and A on the temperature unsteadiness parameter can be shown in this figure.One can demonstrate that the volume fraction and the temperature decline with increasing ε as well as the presence of parameter of unsteadiness A.Figure6illustrates the parameter of magnetic M effect on the velocity, temperature and volume fraction profiles.The temperature only changes near the wall and then it stays uniform through the fluid as it increases.The thickness boundary layer is not under the impact of the magnetic field, due to the volume fraction profile.Therefore, the magnetic field affects the regions that have high concentration of the nanoparticles.The changes near the wall lead to the increment the transfer of heat, which can result in the change in the nanoparticles.Therefore, the temperature changes.Figure7aplots that higher the parameter of the magnetic parameter is, the lower the velocity away from the wall becomes.Physically, we can conclude that when applying the magnetic field inside the boundary layer, it produces a type of resistance force called Lorentz body force that reflects the flow and slows the movement of the fluid.Figure7depicts the impact of Eckert number Ec on the temperature, the velocity, and the heat transfer.Ec statistically affects the fluid temperature because of improving Ec number, the improved viscous dissipation impact on the boundary layer of thermal.The Ec has a strong effect on velocity that we can deduce that the velocity along the surface increases and then becomes steady as η increases.Furthermore, we can see that there is a small change in the boundary layer thickness on the transfer of heat near the wall and then the boundary layer declines to become uniform (no change).Figure8plots the thermophoresis effect Nt and the Eckert number Ec on the friction of skinCf Nu x also, the Sherwood number Sh x .Nt and Ec rise, the skin friction along with Nusselt number decline.In contrast, the Sherwood number rises with rising Nt but declines with the large values of Ec.Figure9displays the impact of Eckert number Ec and Vw on the friction of skin Cf Nu x and Sherwood num- ber Sh x .Cf Nu x rises when rises when Ec and Vw rise, while the Nusselt number and Sherwood number decline with the increasing of Ec and Vw.Figures 10 and 11 display 3D of parameter of magnetic M variation and η, as well, variation of thermophoresis Nt and η on the nondimensional velocity f ′ the surface, respectively.It can be seen in the two cases that the magnetic M and thermophoresis Nt enhance the velocity near the wall.Then, it becomes a constant away through the surface as M, Nt and η increase.Figure12is a plot for the

Figure 7 .
Figure 7. Impact of the variation of the Eckert number parameter Ec for (a) the velocity f ′ , (b) the temperature θ , (c) volume fraction φ and (d) the heat transfer −θ ′ for fixed parameters.

Figure 11 .
Figure 11.Various values of velocity f ′ for varying thermophoresis Nt and η in 3D.

Figure 15 .
Figure 15.Effect of the parameter of mixed convection on Nusselt Number when Scn = 0.22 , Nt = Nb = 0.2 , M = 5 , Pr = 6.2 and Nr = Vw = 0.2 for the variation of unsteadiness parameter A.

Table 2 .
Values of the local Nusselt number concerning various values of Pr, when